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Abstract 

The dynamical response of a relativistic bunch of electrons injected in a pla- 
nar magnetic undulator and interacting with a counterpropagating electro- 
magnetic wave is studied. We demonstrate a resonance condition for which 
the free electron laser (FEL) dynamics is strongly influenced by the presence 
of the external field. It opens up the possibility of control of short wavelength 
FEL emission characteristics by changing the parameters of the microwave 
field without requiring change in the undulator's geometry or configuration. 
Numerical examples, assuming realistic parameter values analogous to those 
of the TTF-FEL, currently under development at DESY, are given for possible 

control of the amplitude or the polarization of the emitted radiation. 
41.60.-m, Radiation by moving charges 41.60.Cr, lasers 
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I. INTRODUCTION. 



Since their first experimental realization JIJ, free electron lasers (FEL) have been one 
of the most promising sources of coherent electromagnetic radiation [@,[|. The physics of 
FEL emission is radically different from that of any other laser sources. In particular, the 
tunability over a broad range of frequencies, and the brightness of its output are difficult to 
achieve in other lasing schemes. On the other hand, its polarization, pulse shape, etc. are 
strongly connected with the geometry of the undulator and hence are inconvenient to modify. 
At the same time, modifications in the typical undulator's physical structure may induce 
new features, e.g., FELs with two magnetic wigglers of different spatial frequencies may 
increase the radiation at higher harmonics suppress the side bands ||, and may allow 
the radiation spectrum || to be controlled. However, a systematic experimental exploration 
of these possibilities is very ackward, if not precluded, due to the difficulty of engineering and 
constructing the modified undulators for every such experiment. It is therefore worthwhile 
to explore theoretically the possibility of achieving control of the amplitude and polarization 
of the emitted radiation, specially at very short wavelengths, without having to alter the 
undulator geometry. 

The basic dynamics of the interaction of free electrons interacting with electromagnetic 
waves has been studied in many circumstances in the past. They range from pioneering 
studies of the radiation of a single electron driven by an electromagnetic wave interac- 
tion of relativistic electrons under general initial conditions with such radiations || , charged 
particle acceleration by simultaneous interaction with an electromagnetic wave and a static 
electric field |10|JTT[] , etc. Another important application of FEL principle is the particle 
acceleration by the inverse mechanism. Particle acceleration by inverse free-electron-laser 
principle has been demonstrated both theoretically ]12| and experimentally [|13| , |14! . In con- 
trast, much less seems to be known for the complementary geometry, in which the electron 
bunch interacts with a counterpropagating electromagnetic wave, perhaps because of the 
absence of acceleration schemes for this case. In this paper we explore the possibility of 
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modifying the electron dynamics in the amplification stage of FEL, by means of the in- 
teraction with a counterpropagating microwave field. It will be shown that, under certain 
conditions, the counterpropagating wave can influence strongly the dynamics of the elec- 
trons inside the undulator. Thus, by a careful choice of the wave parameters control of 
the dynamics can be achieved that could lead to desired FEL radiation properties without 
requiring geometrical changes of the undulators. 



II. ELECTRON DYNAMICS AND PHASE MATCHING CONDITION 

Let us consider a modified FEL configuration as depicted in Fig. [TJ: a free electron is 
injected axially into a linearly polarized magnetic undulator, where an electromagnetic wave 
propagates also axially, in opposite direction, inside a waveguide. The evolution of the elec- 
tron motion, in the combined steady magnetic field of the undulator and the electromagnetic 
wave, is governed by the Newton-Lorentz equation 



d _ 
-p = q 



E + -vx [B u + B 

c 



(1) 



dt 

where B u = f{x)B u sm.(k u x)e z is the undulator's magnetic field. The explicit form of the 
counterpropagating field depends on the waveguide geometry, as well as on the choice of 
a particular transversal mode. Due to the small transversal dimensions of the electron 
bunch used in FEL (of about some tens of microns), only the field at the central axis 
of the waveguide is relevant. The following discussion, therefore, may be applied to any 
waveguide mode of any geometry, provided it has a non vanishing linearly polarized field 
along the central axis, which can be regarded as constant over the whole bunch's section. 
The choice of particular waveguide parameters will influence the quantitative values where 
the interference condition considered below is attained. For concreteness, let us assume a 



TE n o [15] mode propagating in a rectangular waveguide. The explicit forms for the electric 



and magnetic field now read as follows, 



-> k 

E Q = -E g(k w x + uj t) — cos(k c z) sm(k w x + uj t + (f) )e y (2) 

fvr 



B = —E g(k w x + co t) sin(k c z) cos(k w x + co t + (j>o)e x + 

k 

E g(k w x + uj t)-^-cos(k c z) sm(k w x + uj t + 0o)e 2 (3) 

k w = \jk\ — k 2 c being the wave number of the travelling wave, k = ^f, and k c = ™ the cutoff 
wave number of the waveguide (a being the width of the waveguide). f(x) and g{k w x + u^i) 
are considered slowly varying envelopes. 

Without loss of generality (by shifting the time coordinate) we may assume that the 
electron is initially at x — 0, moving along the x-axis with velocity v . Before solving the 
equation of motion numerically we may gain a qualitative insight into the problem by first 
considering the dynamics in a new reference frame in which the electron is initially at rest. 
In the new frame, the undulator's magnetic field becomes a counterpropagating time- varying 
electromagnetic field with 

E'u = -iPfK* + «>'J)B» sin«*' + <t')e y (4) 
B' u = lf(k' u x> + J u t')B u sm(k' u x' + J u t')e z (5) 

where k' u = k u ^, oj' u = k' u v and 7 = l/y/l — (5 2 being the Lorentz factor {(5 = v /c). Note 
that this electromagnetic field has two peculiarities: one, the magnetic and the electric field 
amplitudes do not coincide in their strengths, the electric field being smaller, and two, it 
propagates with a velocity vq = uj' u /k' u < c. A more fruitful way is to reinterpret this field 
as an electromagnetic wave propagating in vacuum with a space dependent phase 



E' u = -^pBJin'y + u'J + (f> u (x')) sin«x' + u'J + <j> u {x'))3 y (6) 

B' u = iB u f{n' u x' + JJ + cf) u (x')) sm(K' u x' + JJ + <p u {x'))e z (7) 

k' u x' + uj'J = k' u x' + u'jf + 4> u {x') (8) 
with < = lv'Jc and <p u {x') =u/ u (±- i) x'. 



On the other hand, the counterpropagating electromagnetic wave in the new reference 
frame becomes 
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E' = -7-E fc ° + cos(k c z')g(k' w x' + u' t') sm(k' w x' + u' t' + (f) )e y (9) 

for 



B' = -E sin(k c z')g(k' w x' + uj' t')cos(k' w x' + uj' t' + o )e x + 

7-E fe °^ + kw cos(k c z')g{k' w x' + u' t') sm(k' w x' + J Q t' + (f) )e z (10) 

with u' Q = j(u> + k w j3c) and k' w = ^{k w + k Q f3). Note that in the strong relativistic case 
(large 7), the effective field observed by the electron can be considered as a TEM wave. 
In addition, the phase velocity of this field approaches c. These two facts, together with 
the condition k c z' ~ 0, which is ensured by the reduced dimensions of the electron bunch, 
permit us to ascribe the effective field acting on the electron in its rest frame to a plane 
wave. 

By inspection of Eqs. (^)-0), one sees that it is possible to derive a phase matching 
condition in which both fields can be seen to have the same frequency in the moving frame, 

K = K 'u } uj kl + k 2 c 



^0 = ^ 



(ii) 



provided that the waveguide has a transverse dimension greater than half the undulator's 
wavelength, a > 4p Note that condition ( |TTD has been calculated for the case of strongly 
relativistic electrons, (3 ~ 1. Although this result is derived for a rectangular waveguide, it 
is worth to stress that this is independent of the particular geometry (which is described 
by the appropiate form of k c ). Note also that the dependence of the frequency of the 
electromagnetic wave on k c permits to attain the same phase matching condition for a 
variety of electromagnetic waves, only by modifying the waveguide geometry. 

As given by Eq. (|TT|) , the phase matching condition is defined only for the temporal 
oscillation. Since the undulator field has a spatial phase dependence, the corresponding 
wavenumber matching will hold only over a certain coherence length such that (t>u{^ C oh) = ^1 

_ X u /2 

Before proceeding further, we may point out that the nature of the electron motion can 
strongly depend on this coherence length and show an interesting disordered behavior when 



the coherence length becomes comparable to the undulator's wavelength. However, here we 
are concerned with the condition in which the coherence length is greater than the total 
undulator length. This condition is easily fulfilled by very high energy electrons. In this 
situation the motion of the electron remains regular. 



III. NUMERICAL INTEGRATION 

As indicated above our objective is to study the possibility of modifying the FEL emis- 
sion characteristics induced by the counterpropagating microwave field, in a configuration 
similar to that being developed in DESY [|16|,[17|. The initial conditions, hence, consist of 
a relativistic electron- bunch entering the undulator in the presence of a very weak seed of 
FEL radiation field, which is assumed to be generated from vacuum noise in a first stage 
of the FEL laser. Since the bunch injection energy is high, the dynamics encloses two very 
different space-time scales, namely, that of the undulator's field and that corresponding to 
the output radiation, which differ typically by a factor j 2 . This disparity becomes a limiting 
difficulty for the numerical integration of the evolution equations, which is usually overcome 
by using the appropriate slowly varying envelope approximations, along with the projections 



on the field cavity modes |18H20|| . In this work we have chosen an alternative procedure |21 



which computes the radiated field from the superposition of the Lienard- Wiechert fields |22| 



emitted from every pseudoparticle (see below) of the bunch. Furthermore, we have preferred 
to integrate the equations in the initial rest frame of the bunch. This allows us to avoid the 
problem associated to the disparity of scales since in this frame the undulator and radiated 
field have similar frequencies. Moreover, in the chosen frame, it becomes readily evident 
that the bunch density is decreased by a factor 7 ^> 1, allowing us to neglect self-fields. 
The large number of electrons per bunch (in our case ~ 10 9 ) in the realistic situation forces 
us to define pseudop&rticles each of which include a few thousands of electrons that are 
assumed to move together. Note that this is the same conceptual philosophy as employed 
in the successful particle- in-cell (PIC) codes for the simulation of plasma dynamics [EBipi. 
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The modulations of the charge density in the system can be modeled either by considering 
the spatially variable distribution of equally charged pseudoparticles, or by a spatially uni- 
formly distributed set of variably charged pseudoparticles. For convenience, we have chosen 
the latter approach in the present investigation. To simulate the velocity and acceleration 
of the pseudoparticles, we have used a relativistic Boris algorithm and to calculate the 



resulting emitted field of the electrons in the forward direction, we have used the well-known 
formula of the far field radiation field amplitude of an accelerated charged particle |22 

q/3 y (t')-{3 x (t')/3 y (t')-l3 y (t)/3 x (t) 



E ra d{t) 



(13) 

t'=t-(R-x(t'))/c 



z (R-x(t>))(l-(3 x (t>))3 

where R is assumed to be large enough. Once the integration is performed, we Lorentz- 
transform the computed quantities to the laboratory reference system. 

IV. COHERENT CONTROL OF FEL RADIATION 

In this section we will theoretically demonstrate the possibility of controlled FEL radia- 
tion through the external electromagnetic wave. The key idea is to consider a counterprop- 
agating wave resonant with the undulator field, in the sense discussed in section [II]. The 
frequency of the wave depends, therefore, on the spatial periodicity of the undulator's mag- 
netic field and on the particular geometry of the waveguide. In our case, we take the 2.73cm 



undulator's wavelength of TTF-FEL at DESY JLy], and a TE\q mode of a rectangular waveg- 



uide of a size of 1.5cm , which is similar to the size of the beam pipe of the FEL at DESY. 
Equation (|ll|) defines the resonant condition for a counterpropagating electromagnetic wave 
in the microwave region with A = 2.99cm when propagating in free space. 

In addition to a resonant frequency, the microwave control of the FEL amplification is 
more effective for the case in which the amplitude of this field in the bunch's rest frame 
equals the amplitude of the electromagnetic wave associated to the undulator's field. At 



present, microwave fields in the GHz range are available with powers up to 100MPF [25 



Although this is already close to the value needed to control optimally the radiation of the 
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TTF-FEL at DESY, we prefer to be conservative and to consider in this paper a tapered 
undulator to reduce the undulator's magnetic field to 25 mT. With this values it should be 
possible to demonstrate the microwave-control experimentally with current technology. On 
the other hand, the state of the art of the microwave generation by the FEL concept allows 
to foresee the availability of brighter sources in the near future [^(J . 

Unless stated otherwise explicitly, the results are based on calculations for an electron 
bunch (of 300 MeV) injected into a 4.5m magnetic undulator, whose characteristics have 
been commented upon above. Our numerical tests show a posteriori that the assumption of 
an initially cold bunch is acceptable. The bunch is described by a spatial sin 2 distribution 
of 20000 particles, 250/im long. Small changes of this number and/or the choice of bunch 
shape is found not to affect the conclusions drawn from the simulations. 

In the following, we consider two cases of microwave control of free electron laser emission. 
First, we will analyze the possibility of suppressing the FEL output by microwave interaction, 
opening ways to control the pulse of the FEL radiation by modulation of the microwave 
amplitude. Second, control of the polarization angle of the FEL radiation by changing the 
microwave polarization. These possibilities are particularly interesting in view of the lack of 
convenient optical elements at very short wavelengths to manipulate these characteristics of 
FEL radiation once they are extracted from the source. 



A. Coherent suppression of radiation 

Let the counterpropagating microwave field be linearly polarized, with the polarization 
vector perpendicular to the direction of the undulator's magnetic field. The undulator field 
and the microwave field may be, then, made to interfere destructively when they have their 
phases properly matched at the a critical value of the amplitude of the microwave field, 

EcrU = W^ Bu (14) 

Note that (like the phase matching condition, Eq. (JIT])) the critical field becomes almost 
independent of the energy of the electron in the highly relativistic case, (3 ~ 1. Thus, a nearly 
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complete destructive interference can occur for the sum magnetic field, or B' T = B' u + B' = 0, 
for all time, in the moving frame. In contrast, because of the asymmetry between the 
magnetic and electric field amplitudes of the (Lorentz-transformed) undulator field, the 
total electric field in the moving frame, on the other hand, does not vanish exactly. A small 
residual electric field E' res = 7(1 — (3)B U remains, which, however, diminishes greatly as the 
electron energy increases. Thus for /3ml, E' res ~ B u /^ 2 , the interference condition for the 
total electric field becomes almost exactly fulfilled. 

Fig. ^| shows the resulting suppression of FEL radiation calculated for different initial 
bunch energies. Note that, as expected, the amplification gain is dramatically reduced as 
the bunch energy increases. This is because the residual electric field , E' res vanishes with 
increasing energy of the bunch. Note that, for the higher energies, the microwave field 
reduces the gain by nearly 3 orders of magnitude. 

For the destructive interference mechanism to be effective in practice, it is required a 
constant tt phase-difference between the undulator and the microwave fields, as seen in the 
electron bunch's reference system. Note that ensuring an initial ir dephase requires a certain 
control of the bunch's conditions before injection, since the bunch must enter the undulator 
when the microwave phase is opposite to the undulator's. The required constancy is ensured 
by (|12D since the coherence length is greater than the undulator's dimension for the range 
of bunch energies assumed here (e.g., TTF-FEL). 

To analyze the sensitivity of the coherent suppression effect 

against fluctuations in the initial field dephasing, we have performed a series of calcu- 
lations in which the bunch's initial position against the undulator vertex is changed. The 
shift in the initial position is directly related to the time delay of the bunch to reach the 
undulator and, therefore, to the initial dephase between the undulator and microwave fields 
(in the rest frame of the bunch). The results of calculations with a 300MeV, 250/im bunch 
are presented in Fig. |3] which shows the amplification vs. fluctuations in the bunch position. 
It can be seen that the gain suppression effect is robust against fluctuations less than 1mm 
that is well above the usual experimental uncertainty. 
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B. Control of polarization of FEL radiation 

Polarization control is of particular interest for very short wave FEL radiation. This 
can be achieved in the same configuration by rotating the polarization of the microwave by 
a certain angle from the undulator's plane of polarization. From the theoretical analysis 
above we expected that, in general, the emitted radiation will be elliptically polarized. 
The ellipticity would depend on the initial dephase of the fields as well as their relative 
amplitudes. We show in Fig. |]a the calculated change of ellipticity of FEL output vs. the 
relative angle between the polarization plane of the undulator and that of the microwave 
(chosen to be linearly polarized). The amplitudes of the undulator and microwave fields in 
the rest frame of the bunch are chosen to be comparable while the initial dephase is set to 
7r/2. The resulting ellipticity of the emitted radiation is found (Fig. |]a) to change from the 
linear to the circular polarization. The results of the various cases presented in this figure are 
summarized in Fig. |]b in terms of the tilt-angle of the major axis of the polarization-ellipse 
of the emitted radiation. 

V. CONCLUSION. 

Coherent modifications of FEL radiation induced by a counterpropagating electromag- 
netic wave interacting with an electron bunch in a magnetic undulator are studied. A 
phase-matching condition between the undulator field and an external microwave field in 
the rest frame of a relativistic electron bunch is derived. This condition is found to only 
depend on the geometry of the problem. It is found that possible control of both the ampli- 
tude and the polarization of the FEL radiation (including very short wavelengths) could be 
achieved without having to alter the undulator's geometry, by simply varying the incident 
microwave field. Results of concrete numerical simulations assuming realistic FEL parame- 
ters (corresponding to that of TTF-FEL, currently under development at DESY) are given, 
and their robustness against small fluctuations in initial conditions is illustrated. 
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FIGURES 

FIG. 1. Schematic diagram of the modified FEL amplifier configuration used throughout this 
paper. An electron-bunch is injected into the linearly polarized magnetic undulator in the presence 
of a counterpropagating electromagnetic wave. 

FIG. 2. Dependence of the FEL amplification factor versus the initial bunch energy, in the 
presence of the microwave field (dashed line) and in its absence (solid line). 

FIG. 3. Dependence of the amplification factor versus the initial coordinate of the bunch; x = 
corresponds to the initial position for which the undulator and the microwave fields in the rest frame 
of the bunch have opposite phases. 

FIG. 4. (a) Ellipticity of the amplified FEL radiation versus the relative angle between the 
planes of polarization of the microwave and the undulator fields, (b) Tilt-angle of the major axis of 
the polarization ellipse of the FEL radiation with respect to the polarization plane of the undulator 
field, for the same cases as considered in (a). 
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